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ABSTRACT 

We have calculated the 0(a ew Mj- /M^) supersymmetric electroweak corrections to 
single top quark production via qq' — > tb at the Fermilab Tevatron in the minimal 
supersymmetric model. The supersymmetric electroweak corrections to the cross 
section are a few percent for tan (3 > 1, and can exceed 10% for tan/3 < 1. The 
combined effects of SUSY electroweak corrections and the Yukawa corrections can 
exceed 10% for favorable parameter values, which might be observable at a high- 
luminosity Tevatron. 
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1. Introduction 

The top quark has now been discovered by the CDF and DO collaborations at the Fer- 
milab Tevatron[l]. Measurements of its mass are 176 ± 9GeV and 170 ± 18 from CDF and 
DO, respectively, and the world- average value of the top quark mass from Run I at the 
Tevatron was recently reported to be 175 ± QGeV, based on roughly 100 pb~ x of data[2]. 
At the Tevatron the dominant production mechanism of the top quark is the QCD pair 
production process qq — > tt[3\. However, because the top quark is so heavy, electroweak 
production of single top quarks can become significant, particularly at the next Tevatron 
Run-II. With y^s = 2 TeV and an integrated luminosity of 2fb~ l one can expect, in the 
Standard Model(SM), that for a 175GeV top quark, there will be about 1.4 x 10 4 tt pairs 
and 5 x 10 3 single top events produced[4], which is about 35% of the total tt rate. After 
taking into account the 6-tagging efficiency and the detection efficiency [5], there are about 
1000 single-6-tagged tt pairs in the Z+jets sample, 100 in the dilepton sample, and 250 single 
top events in the Z+jets sample available for testing various properties of the top quark. 
Even with fewer events, single top production processes are important because they involve 
the electroweak interaction and, therefore, can probe the electroweak sector of the theory, 
in contrast to the QCD pair production mechanism, and provide a consistency check on the 
measured parameters of the top quark in the QCD pair production. At the Tevatron sin- 
gle top quarks are produced primarily via the lU-gluon fusion process [6] and the Drell-Yan 
type single top process, qq' — > W* — > tb(W* process) [7], which can reliably be predicted in 
the SM. The theoretical uncertainty in the cross section is only about a few percent due 
to QCD corections[8]. As analysed in Ref.[9], the statistical error in the measured cross 
section for the W* process at the Tevatron will be about ±30%; however, a high-luminosity 
Tevatron would allow a measurement of the cross section with a statistical uncertainty of 
about 6% [9]. At this level of experimental accuracy a calculation of the radiative corrections 
in the SM is necessary and effects beyond the SM, for example SUSY corrections, should 
also be considered. 

In Ref.[9] the QCD and Yukawa corrections to the W* process have been calculated 
in the SM. In a previous paper [10] we calculated the Yukawa corrections to this process 
from the Higgs sector in the general Two-Higgs-Doublet Model(2HDM) and the Minimal 
Super symmetric Model(MSSM)[ll] and found that the corrections can amount to more than 
a 15% reduction in the production cross section relative to the tree leve result in the 2HDM, 
and a 10% enhancement in the MSSM. However, in the MSSM, in addition to these Yukawa 
corrections from the Higgs sector, the Supersymmetric (SUSY) corrections due to super 
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particles (sparticles) should also be taken into account. The dominant virtual effects of 
sparticles arise from the SUSY electroweak corrections of order 0{a ew Mf / M^) and the 
SUSY QCD corrections of order 0(a s ) which arise from loops of charginos, neutralinos and 
squarks, and gluinos and squarks. It is well-known[12] that the anomalous magnetic moment 
for a spin 1/2 fermion vanishes in the SUSY limit and away from the SUSY limit there is a 
partial cancellation. Therefore, in general, one can expect the Yukawa corrections from the 
Higgs sector and the SUSY electroweak corrections from virtual charginos and neutralinos 
to cancel to some extent. 

In this paper we present the calculation of the 0(a ew M 2 /Myy) SUSY electroweak cor- 
rections at the Fermilab Tevatron in the MSSM, and show the combined effect of including 
all the 0(a ew Mf / M$y) terms from both Yukawa corrections and SUSY electroweak correc- 
tions. The 0(a s ) SUSY QCD corrections will be given elsewhere [13]. The paper is organized 
as follows: In Sec. II we present the analytic results in terms of the well-known standard 
notation of one-loop Feynman integrals. In Sec. Ill we give some numerical examples and 
discuss the implications of our reults. 



2. Calculations 



2.1 Formalism 

In our calculations we used dimensional regularization to control all the ultraviolet di- 
vergences in the virtual loop corrections and we adopted the on-mass-shell renormalization 
scheme[14]. Including the 0(a ew Mf / M^) electroweak corrections, the renormalized ampli- 
tude for qq' — > tb can be written as 

M ren = M + SM SUSY - EW (1) 

where M is the tree-level matrix element and §m susy ~ ew represents the SUSY electroweak 
corrections. The tree-level Feynman diagram for single top quark production via qq' — > tb is 
shown in Fig. 1(a). The amplitude M is given by 

q 2 1 

M ° = 1 2 s - M^ &h^PMptMp^PLvipi). (2) 
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The amplitude §M SUSY - EW is given in the following section. Here p\ and p 2 denote the 
momentum of the incoming quarks q and q', while p^ and p 4 are used for the outgoing t and 
b quarks, and s is the center-of-mass energy of the subprocess. 
The renormalized differential cross section for the subprocess is 

da ^ — M?^; 12 
a cos 9 32tts 2 ^ 

where # is the angle between the top quark and incoming quark. Integrating this differential 
cross section over cos 6 one obtains the cross section for subprocess 

a = a + Act (4) 

where the tree-level cross section ct is given by 



and Act represents the SUSY electroweak corrections. 

The total hadronic cross section for the single production of top quarks via qq 1 can be 
written in the form 

a(s) = Y,j dx 1 dx 2 a lJ (x 1 x 2 s,M^fi 2 )[f l A (x 1 ,fi)ff(x 2 , f i) + (A <- B)], (6) 



where 



s = (P1 + P2) 2 , (7) 
s = xix 2 s, (8) 
Pi = X1P1, (9) 



and 



p 2 = x 2 P 2 . (10) 

Here A and B denote the incident hadrons and Pi and P 2 are their four-momenta, while 
i,j are the initial partons and x\ and x 2 are their longitudinal momentum fractions. The 
functions ff and f? are the usual parton distributions[15,16]. Finally, introducing the 
convenient variable r = X\x 2 and changing independent variables, the total cross section 
becomes 

i j j tq t s ar 
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where r = (M t + M b ) 2 /s. The quantity dL^/dr is the parton luminosity, which is defined 
to be 

dJ ^= f —\f?(xi,li)fj i (T/x 1 ,n) + (A~ B)} (12) 
dr Jt x\ j 



2.2 SUSY electroweak corrections 

The SUSY electroweak corrections of order a ew M 2 / 'M^ to the process qq' — > tb arise 
from the Feynman diagrams shown in Fig.l(b)-(g). The matrix element for these corrections 
can be written as 



5M SUSY-EW = 5M SUSY-EW + + ^ 



box 



(13) 



where SM^^ Y ~ EW represents corrections arising from the self-energy diagrams and vertex 
diagrams [Fig.l(b)-(e)], while 5Mb ox and SM£ ox correspond to the box diagram [Figl.(f)] and 
crossed box diagram[Figl.(g)], respectively. 5M^/£ Y ~ EW is given by 



SUSY-EW 
Wtb 



. 9 Z 



2 



v(P 2 )^P L u(P 1 )u(P 3 



^Pl^SZI + l -5Z b L + Et) 



+PgP L E% + P£P L E E v(P 4 ) 



(14) 



The renormalization constants and form factors in Eq.(14) are 



5Z\ 



8Z b L 



E L 2 



and 



8tt 2 



+M^(|^o| 2 +|L^| 2 )C7^ ) 



_9 \2r2| T/ \2 M 2 r ^i b i) 



E? = 



\ t V; 2 R- ti Ol] \L^oM^c 12 + R m M t (c 12 + c 22 ) 



(15) 
(16) 



^- 2 ^V* 2 R k (i^oOg M, ; M,oc + L^Og M, ; M,(c + c 12 ) 

+^00^ M^oM t c 12 + Ri^pl]\Ml{c 22 - c 23 ) + s(c 12 + c 23 ) + 2c 24 - -]} ,(17) 



E 



3 L = ^^^ V h R u l L nxoO«*M- t (c + c u ) - R-^pi; M t (c 12 + c 23 ) 



(19) 
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where sums over k are implied and the functions Q,(P 4 , P 3 , M-+, Mf., M^o) are the three- 
point Feynman integrals [17]. The functions F^{ k \ and constants in the above equations 
are defined as 



rp(ijk) 
± n 

r<{ijk) 

Xt 



dyy n log 



dy 



™%y(y - !) + - y) + m ^ 



fJL- 

n+1( i-y) 



o m 2 ty(y - I) + m 2 Al - y) + mly' 



M w sin /5 ' 
cos On, L n 



Ra, — 



^IX? 
^?2X° 



— sin 9q, 

cos 



and 



sin ^ R , Q2 — ' ! f ' f 7 " 

Aj cos — Cj sin ^, 
— A,- sin 9q — Cj cos Oq, 
—A* sin + Bj cos 0,, 
—A* cos 9q — Bj sin 0^, 



0% = ^n; 3 u j2 + N; 2 u n , 



where 9 q is the mixing angle of squark q, and 



A, = 



C 3 = 



gm q N, 



J"4 



s,- = c* + 



gN' 



J 2 



2m w sin/3' J ~ 3 ' 2CW 



2 gSy^^ lJf 



- e N * - 
3 Jl 3 CV 

NjiCw + Nj 2 Sw, 



J2' 



Ar; 2 = -N^Sw + NjtCw. 



(20) 
(21) 

(22) 

(23) 
(24) 
(25) 
(26) 
(27) 
(28) 

(29) 
(30) 

(31) 

(32) 

(33) 
(34) 



In the above Sw = sin 9w and Cw = cos0vi/. The chargino masses M,- and matrix elements 
Vij depend on parameters M,/i and tan/5, whose expressions can be found in Ref.[ll]. The 
neutralino masses M^o and matrix elements are obtained by diagonalising the matrix 
Y [11]. Given the values of the parameters M,M',fi and tan/3, the matrix iV and M^o 
can be obtained numerically. Here, \x is the coefficient of the Hi — H 2 mixing term in the 
superpotential and M and M' are the masses of gauginos corresponding to SU(2) and U(l), 
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respectively. With the grand unification assumption, i.e. SU(2) x U(l) is embedded in some 
grand unified theory, we have the additional relation M' = I^M. 
The box diagram amplitude 5M box is 



SM, 



box 



iu(P 3 ) [(f*P R + f 2 P L + ft f 4 P R + fl f 4 P L ) u(P 1 ) v(P 2 ) f z P L 
+ (ft-fPn + fh'Pi) «(Pi) v{P 2 )^P L ] v(P 4 ) 



(35) 



Here the form factors f 435 are 

,2 



fi 



ft 
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L m M t (D 12 - D 13 + D 22 - D 26 ) 
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8 V27r 2 XtVjlV V Rii L ~ qi L *«ix°k L *«° ( Dl2 ~ Dl3 + ° 24 ~ J ° 25 )' 



ft - - ^tVjiV^Rf. Lq t L*~ -o L^o D 27 

and /| j4i6 can be obtained through the permutation 

fb rb 

i2,4,6 ~~ J 1,3,5 



(36) 
(37) 

(38) 



Hx° k t;x°' qix ° qi x° 



(39) 



The sums over i, j, k, I are implied and the functions Dij(P 4 , P 3 , —Pi, M~ to M^,o, M^) are 
the four-point Feynman integrals [17]. 

The amplitude for the crossed box diagram 5M£ ox is 

5M c box = -iu(P 3 ) [A C P L + f c 2 P R + / 3 C APl + / 4 C APr] «(P 2 ) v(Pi)Plv(P4), (40) 
where the form factors are 

„2 



"i^^ W^^^x2 M x; i M ^x° fc (A) + A 2 ) + M x oL^„P 



and 



ft = -^ At ^^ L ^« M x 3 +i '-^^ + ^ 



(41) 
(42) 
(43) 

(44) 
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The sums over k, I are again implied and the functions D^P^, P 3 , —P2, Mf., M^o, M^] 
are the four-point Feynman integrals [17]. The constants L'- ^ ,R'~ ^ are defined by 



and 



with 



and 



L kx° = ^cos^-Cjsin^, (45) 

L' h9 n = -^sin^-Cjcos^-, (46) 

R'^o = -A'* sin 0,- + £;.cos0,-, (47) 

^kyo = -4* cos h ~ B 'i sin (48) 



J 2m w cos/3' J J 2CV' 1 J 



05 = 4^ + 1^. (50) 



3. Numerical results and conclusion 

In the following we present numerical results for the corrections to the total cross section 
for single top quark production via qq' — > tb at the Fermilab Tevatron with ^/s = 2 TeV. In 
our numerical calculations we used the MRSG parton distribution functions[16] and chose the 
scale [i = V§. Also we neglected SUSY corrections to the parton distribution functions. For 
the parameters involved, we chose M z = 91.188GeV, M w = 80.33GeV,M t = 175GeV, M b = 
5GeV and a ew = 1/128. Other parameters were determined as follows: 

(i) The upper bound on tan/3; viz, tan/3 < 0.52GeV~ 1 M H +, was determined from data 
on B — > rz/X[18]. The lower limits on tan/3 are tan/3 > 0.6 from perturbative bounds [19] 
and tan/3 > 0.25 (for M t = 175GeV) from perturbative unitarity[19]. We limited the value 
of tan /3 to be in the range of 0.25 to 5, as larger values of tan /3 are not interesting, although 
allowed by the current data[18], since the effects are negligibly small. 

(ii) For the parameters M^ R , M^ L , tan /3 and M LR = A t +(j, cot (5 in top squark (stop) mass 
matrix [20] 

2 _( Mf L +m 2 t + 0.35 cos(2/3)Mf -m t (A t + \l cot (3) 

1 ~ 1 L -m t (A t + /I cot (5) M~ +m 2 t + 0.16 cos(2/3)M| 
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we assumed Mf = M^ L . There are then three free parameters in the stop sector and we chose 
the mass of the lighter stop , M^ R and tan (3 to be the three independent parameters. 
The best current lower bound on the stop mass is 55GeV coming from LEP, operating at 
s/s = 130 — 140 GeV[21]. We conservatively took the lower bound to be 50 GeV for m^. For 
the other squarks; i.e., q,q' and b, we neglected the mixing between left- and right-handed 
states and assumed = Mg 2 = = = = which was then determined 
by[20] 

ml =ml + M? L + cos(2/?)(-i + \s 2 w )M 2 z (52) 

(hi) For the parameters M, M',/i and tan/5 in the chargino and neutralino matrix, we 
put M = 200 GeV, ^ = -100 GeV and then used the relation M' = f^M [11] to determine 
M'. 

Some typical numerical calculations of the SUSY electroweak corrections are given in 
Figs.2-4. 

Fig. 2 shows the SUSY electroweak correction Aa/a as a function of lighter stop mass 
Mf l5 assuming tan (3 = 1 and M LR = m t . The correction is only a few percent of the tree- 
level value a and is quite sensitive to the lighter stop mass. There are two peaks at about 
M h = 75 GeV and 67GeV due to the fact that m t = 175 GeV, M^o = (100, 107, 128,221) 
GeV and the threshold for open top decay into a neutralino and a lighter stop is crossed in 
these regions. 

Fig. 3 gives the SUSY electroweak correction as a function of M LR , assuming tan (3 = 1 
and Mi x = 60 GeV. This correction is also sensitive to M LR . With increasing M LR the mass 
splitting between the two stop quarks increases. Since we fixed the mass of ti, the mass 
of t 2 then increases with M LR . Furthermore, with an increase of M LR) M- tL increases and 
thus the sbottom masses also increase, as seen from Eq.(52). Since we assumed the masses 
of the squarks q~i,q~2,Qi an d q' 2 are degenerate with the sbottoms, these masses then also 
increase with M LR . So, with an increase of M LR) all the squark masses except the lighter 
stop increase and their virtual effects decrease due to decoupling effects. From Fig. 3 one sees 
that for M LR > 200 GeV the magnitude of the correction drops below one percent. 

In Fig. 4 we present both the SUSY electroweak correction and the Yukawa correction [10] 
as a function of tan (3. For the SUSY electroweak correction we assumed M LR = m t . Since 
both these corrections are proportional to m2 s - n ^ , they can be very large for small tan/3. 
From Fig. 4 one sees that the SUSY electroweak correction exceeds -10% for tan/3 < 0.5, 
and tends to infinity when tan/5 tends to zero. As in the case of the Yukawa corrections, 
the SUSY electroweak corrections are negligibly small for tan/5 > I. Also, comparing the 
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SUSY electroweak correction with the Yukawa correction in Fig.4, one notes that the SUSY 
electroweak correction and Yukawa correction have opposite signs, and thus cancel to some 
extent. If the lighter stop quark has the same mass as the charged Higgs boson, the cancel- 
lation is appreciable. However, as seen in Fig.4, if the charged Higgs is much heavier than 
the lighter stop quark, the magnitude of Yukawa correction is much smaller than the SUSY 
electroweak correction and there is very little cancellation. In such a case the combined 
effects can exceed -10% for tan/3 < 1. 

To summarize, the combined effects of SUSY electroweak corrections and the Yukawa 
corrections can exceed 10% for favorable values of the parameters. Since the cross section 
for single top production can be reliably predicted in the SM [9] and the statistical error 
in the measurement of the cross section will be about 6% at a high-luminosity Tevatron[9], 
these corrections may be observable; at the least, interesting new constraints on these models 
can be established. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant No. DE-FG02-91-ER4086. 
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Figure Captions 



Fig.l Feynman diagrams for the SUSY electroweak (SUSY EW) corrections: (a) tree- 
level, (b)-(h) one-loop corrections. 

Fig. 2 The SUSY electroweak (SUSY EW) correction Aa/a as a function of M^, assum- 
ing tan (3=1 and = m t . 

Fig. 3 The SUSY electroweak (SUSY EW) correction as a function of M LR , assuming 
tan f3 = 1 and M~ tl = 60 GeV. 

Fig. 4 The SUSY electroweak (SUSY EW) correction and the Yukawa correction as a 
function of tan (3. M LR = m t was assumed for the SUSY EW correction. 
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